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Abstract 

By using the method of density-matrix renormalization-group to solve the different spin-spin corre- 
lation functions, the nearest-neighbouring entanglement(NNE) and next-nearest-neighbouring entangle- 
ment(NNNE) of one-dimensional alternating Heisenberg XY spin chain is investigated in the presence of 
alternating nearest neighbour interactions of exchange couplings, external magnetic fields and next-nearest 
neighbouring interactions. For dimerized ferromagnetic spin chain, NNNE appears only above the critical 
dimerized interaction, meanwhile, the dimerized interaction effects quantum phase transition point and im- 
proves NNNE to a large value. We also study the effect of ferromagnetic or antiferromagnetic next-nearest 
neighboring (NNN) interactions on the dynamics of NNE and NNNE. The ferromagnetic NNN interaction 
increases and shrinks NNE below and above critical frustrated interaction respectively, while the antifer- 
romagnetic NNN interaction always decreases NNE. The antiferromagnetic NNN interaction results to a 
larger value of NNNE in comparison to the case when the NNN interaction is ferromagnetic. 
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I. INTRODUCTION 



Entanglement is one of the most profound features of quantum mechanics and has been 
considered as an important resource in quantum information processing (QIP) including 
teleportationj^l cryptography and secure direct communication, and quantum communication 
and computation. In order to realize quantum information process, great effort has been devoted 
to the generation of entanglement in linear optics,'^] cavity QED t^"^' and ion trap'^'^'^^ schemes. In 
recent years much attention has been focused on the spin system with Heisenberg exchange inter- 
action, which is a typical quantum system, and the entanglement properties have been extensively 
investigated in spin systems, such as Ising model'^^l and isotropic and anisotropic Heisenberg 
models. However; as far as we know, most discussions mentioned above merely focused on 
the models with nearest neighbour interaction of exchange couplings, while next-nearest neigh- 
bour interaction has not been taken into account. Dimerized systems and frustrated systems with 
next-nearest- neighbour (NNN) interactions play an important role in condensed matter the- 
ory, and entanglement in dimerized and frustrated systems has been considered. Sun et al'^^' 
investigate entanglement properties in dimerized and frustrated spin-one models by applying the 
concept of negativity. In Ref.[20], entanglement is studied in an open alternating chain of nuclear 
spins s = 1/2 with spin-spin couplings in an external magnetic field under the thermodynamic 
equilibrium conditions. Gu et al'^^1 have investigated entanglement in frustrated spin-half Heisen- 
berg chains. Chen Yan et al'^^1 analyzed sublattice entanglement and quantum phase transitions in 
antiferromagnetic spin chains, Chen Shu et al'^^^ study the fidelity and quantum phase transition for 
the Heisenberg chain with next-nearest-neighbouring interaction. In fact, some one-dimensional 
and two-dimensional compounds (CuGeOs, NaV20^) have manifested such interactions. There- 
fore, it is worthwhile to include next-nearest neighbouring interaction in the studies of spin chain 
entanglement. 

The ID Heisenberg model is a simple but realistic and extensively studied in solid state system. 
Recently, Osterloh et al '^^^ examined the entanglement between two spins of position i and j in 
the spin chains for the pure case, Huang et al,^^' Osenda etal'^^^ and we^^^' have demonstrated that 
entanglement can be controlled by introducing impurities into the systems. Meanwhile, the entan- 
glement shares many features with quantum phase transition (QPT) for a many -body system. 
QPT, which occurs at absolute zero temperature and is purely driven by quantum fluctuation, is 
the structural change in the properties of the ground state. The associated level crossings lead 
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to the presence of non-analyticities in the energy spectrum. Therefore, the knowledge about the 
ground-state entanglement, the nonlocal correlation in quantum systems, is considered as a key 
to further understand QPT. In this paper, we study the pairwise ground state entanglement be- 
tween the nearest-neighbor sites and that of the next-nearest neighbouring sites in one-dimensional 
s — ^ Heisenberg XY spin chain with dimerised exchange couplings and next-nearest neighbour 
coupling, to our knowledge, which has not been reported before. The present study will help us 
to further understand the behavior of the entanglement and QPT in one-dimensional alternating 
Heisenberg XY model with next-nearest neighbour interactions, and this model can display a va- 
riety of interesting physical phenomena. More interestingly, we can control or manipulate the 
NNE and NNNE in QPT point with the help of dimerised exchange couplings, dimerised external 
magnetic fields and frustrated exchange interaction. 



II. SOLUTION OF THE XY MODEL AND SPIN-SPIN CORRELATION FUNCTIONS 

We consider one-dimensional alternating Heisenberg XY model of N spin-| particles with 
nearest-neighbour interactions and next-nearest neighbour interactions. One-dimensional Hamil- 
tonian can be written as 

AT 1 _ ^ ^ 

^ = - ^ E -^M+i^+i - ^ E '^v+i'^r'^r+i - E ^^^^ 

i=l i=l i=l 

where Ji^i+i and are alternating exchange interaction and nearest-neighbour interactions re- 
spectively, hi is the strength of external magnetic field on site i, a^'V'^ are the Pauli matrices, 7 is a 
dimensionless parameter characterizing the anisotropy of the model and N is the total number of 
sites. Furthermore, the periodic boundary conditions satisfy cr^+i = erf, cr^+i = erf, cr^+i = cr^. 

The spectrum of this Hamiltonian can be determined exactly by a straightforward applica- 
tion of the standard methods. The first step in the procedure is to perform a Jordan-Wigner 
transformationt'^^l by introducing Fermi operators Cj, Cj, as a result, the Hamiltonian (1) is mapped 
in the free fermion Hamiltonian 

N ^ 1 

^ = - E J^,^+l[i4c^+l + 7C+C+ J + h.c] - 2 E hi{cic, - -) 
i=l i=l 
N 

- E '^i,i+2[i4ci+2 + 144+2) + h.c] (2) 



i=l 



In this paper, the exchange interaction and external magnetic field have the form J2i-i,2i — Ji — J, 
J2i,2i+i = ^2 = aJ, h2i,2i+i = h2i-i,2i{h) = i = 1, 2, • • • , y) respectively, where a = J2/ Ji 
introduces the dimerised parameter, the ratio of next-nearest (Ji,i+2 = ^3) to nearest-neighbor 
interaction coefficients k = J3/ Ji is called the frustration parameter. J < and J > correspond 
to the antiferromagnetic and the ferromagnetic cases, respectively. For J3 = 0, Eq. (2) reduces 
to one-dimensional altemating(dimerised) Heisenberg model, for Ji = J2, it reduces to one- 
dimensional frustrated Heisenberg model. 

The density-matrix renormalization-group (DMRG)method[^^'^^lis applied to obtain the spin 
correlation function. The two-point reduced density matrix obtained by tracing the full density 
matrix of the system over all sites except the pair {i, j} has the form 

k=x,y,z 

In this part, we choose the concurrence defined by Wootters '^^1 as a measurement of the pairwise 
entanglement. For a pure or mixed state of two qubits described by the density matrix p, the 
concurrence C may be calculated explicitly as 

C{p) = max(0, Ai - A2 - A3 - A4) (4) 

where Ai, A2, A3, A4 are the eigenvalues in descending order of the spin-flipped density operator 
R, which is defined by R = ^J~^Jpf\fp , where p = {ay ® ay)p*{ay ® Oy), p denotes the complex 
conjugate of p, ay is the usual Pauli matrix. All the matrix elements in the density matrix can be 
calculated from the different spin-spin correlation functions. 

III. RESULTS AND DISCUSSIONS 

In this paper, we focus our discussion on the transverse Ising model with 7 = 1 and introduce 
the dimensionless parameter A — J/2h. The goal of our present study is to find dynamic 
characteristics of the entanglement in varying the dimerised parameter of exchange couplings, 
external magnetic fields and NNN interaction. Figure 1 displays the diagram of the Ji — J2 — J3 
model which can be considered either as a linear chain(a) or a zigzag chain(b). The alternating 
nearest-neighbouring interaction is Ji and J2 respectively, while the next-nearest-neighbouring 
interaction is parameterized by J3. We consider three different physical models, i.e, dimerized 
ferromagnetic spin chain, frustrated ferromagnetic-ferromagnetic(F-F) spin chain, frustrated 

4 



ferromagnetic-antiferromagnetic(F-AF) spin chain. First, we examine the dimerized Heisenberg 
chain, this model is characterized by an alternation of strong and weak bonds between two 
nearest-neighbouring spins. Figure 2 shows the nearest neighbouring concurrence C(l,2) 
and C(2,3), which correspond to two nearest-neighbouring spins coupled by bonds Ji and J2 
respectively, as a function of the reduced coupling constant A at different values of the dimerised 
parameter a of exchange couplings, parameter P of external magnetic fields with the system size 
N =59. Numerical results in Figs.2(a) and 2(c) show that the concurrence C(l,2) decreases when 
the dimerised parameter a of exchange couplings increases, while for the concurrence C(2,3) the 
situation is opposite, and the pair of qubits (1,2) and (2,3) have the same pairwise entanglement 
for the homogeneous chain a = 1. In fact, in the case of weak alternation of the nearest- neighbor 
exchange coupling, we have here the dimerised spin chain which can be considered qualitatively 
as a set of non-interacting spin pairs. In the limit of strong dimerization, the concurrence reaches 
a maximum value. Moreover, spin 2 can be entangled both with spin 1 and with spin 3. Since 
spin 2 is strongly entangled with spin 1, the entanglement of spins 2 and 3 is weaker. Obviously, 
the increase of exchange interaction J2 suppresses the entanglement of spins coupled by Ji. The 
increase of J2 can enhance the entanglement between spins coupled by J2. This explains the 
behavior of the concurrence displayed in Fig. 2. The effect of the alternation of external magnetic 
field P is also shown in Figs.2(b) and 2(d). However, different from the effect of the exchange 
couplings, the concurrence (1,2) and (2,3) have similar behavior and increase firstly and then 
decrease with increasing the value of the parameter /3, suggesting a critical alternation of external 
magnetic field Pc, where the maximum concurrence occurs, must be exist. 

From Fig.2, we can see that the alternate interaction plays an important role in enhancing the 
nearest neighbouring concurrence. The maximum value of the concurrence between neighboring 
sites does not occur at the critical point. The reasons are based on the properties of shared entan- 
glement to expect that this maximum should occur away from the critical point. Entanglement 
sharing is relevant to the quantum phase transition in the transverse Ising model as it provides 
a fundamental bound on the amount of entanglement that may be distributed among the other 
sites, which means that as the overall entanglement in the lattice is increased, some sites become 
disentangled. However, as we all know, the maximum value of the concurrence between the 
next-nearest-neighbouring sites does occur at the critical point A = 1, so it is necessary to discuss 
in this subsection the effect on the next-nearest-neighbouring entanglement(NNNE). Figure 3(a) 
and 3(b) show the change of concurrence C(l,3) as a function of A for different values of alternate 
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interactions a and /3 with k = 0, i.e. in the absence of NNN interaction. We can see that there is 
no entanglement at a = 0.5 and /3 = 0.5. The weak alternate interactions a and /3 suppress and 
finally completely destroy the NNN concurrence. The dashed line in Fig.3(a) depicts at the critical 
point, A = 1, of the transverse Ising model(a = 1), which maximum next-nearest-neighbouring 
entanglement occurs, there is a fundamental transition in the structure of the ground state. This is 
consistent with the former result. As a increases the NNN concurrence tends to increase faster, the 
peak value induced by the NN interaction increases and the A^, where concurrence approaches 
a maximum, shifts to left very rapidly. The dimerized interaction does assist the entanglement 
of formation because of the inside alternation. Therefore, we can further understand the relation 
between the entanglement and quantum transition. In Figs. 3(c) and 3(d), we plotted our numerical 
results for the threshold NN alternate interactions a and /3 for the case of next-nearest-neighboring 
concurrence. Only above the critical a and NNN concurrence appears. 

For the frustrated model, there is a competition between the NN and NNN interactions. So 
the model displays some special features of entanglement. In Fig.4, we give the results of the 
NN and NNN concurrence as a function of the parameter A for different frustrated interaction k. 
The NNN interaction has a different frustration effect on C(l,2). For frustrated ferromagnetic- 
ferromagnetic(F-F) spin chain /t > 0, by the comparison among the different curves in Fig.4(a), it 
is interesting to find that the entanglement peak between the nearest neighbours in the dotted line 
increase to a value larger than those in the solid line. With the increasing of k, in the dash-dotted 
line, the concurrence decreases. We can imagine there must be a critical frustration strength (kc), 
below Kc, NNN interaction enhances entanglement, while above Kc, NNN interaction shrinks 
entanglement. This behaviour is due to the energy level cross at the point k = 0.5, as seen clearly 
from Fig.4(a). Before and after the critical point, the entanglements display distinct behaviours. 
However, for frustrated ferromagnetic-antiferromagnetic(F-AF) spin chain k < 0, the situation 
becomes different due to the frustration, with the increase in absolute value of k and weak 
parameter A, antiferromagnetic NNN interaction always shrinks entanglement, which is different 
from the results obtained from Fig.4(a). Another important character revealed by Fig.4(b) is 
that, as parameter A increases, the concurrence increases and tends to move to infinity(> 0) by 
increasing the value of the parameter k, the strong k is helpful to keeping the better entanglement. 
It is the frustrated interaction that leads to considerable different evolutions of the entanglement, 
hence the entanglement is rather sensitive to any change with the NNN interaction. From Figs. 4(c) 
and 4(d), we can see that the height of the peak increases with the increase of antiferromagnetic 
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K, which is a characteristic feature of one dimensional Heisenberg antiferromagnetic chain. This 
suggests that the antiferromagnetic component becomes significant for larger k,. Meanwhile, the 
antiferromagnetic NNN interaction results to a larger value of NNNE in comparison to the case 
when the NNN interaction is ferromagnetic. In order to explain these results it is necessary to 
take into account the influence of the competing roles played by NN and NNN interaction on the 
entangled states. Thus, by adjusting NNN interaction constant, one can control or manipulate the 
NN and NNN entanglement. 

In summary, from the above analysis, it is clearly noted that dimerised exchange couplings, 
dimerised external magnetic fields and frustrated exchange interaction, which play the competing 
roles in enhancing quantum entanglement, have a notable influence on NN and NNN concurrence 
in one-dimensional s = ^ XY spin system. The NN and NNN concurrence exhibits some 
interesting phenomena. For dimerized ferromagnetic spin chain, NNNE appears only above the 
critical dimerized interaction, which effects quantum phase transition point and increases NNNE. 
For a ferromagnetic-antiferromagnetic frustrated spin chain, the NNN interaction is predominant 
in the competing role and can enhance the entanglement to a steady value. The ferromagnetic 
NNN interaction increases and shrinks NNE below and above critical frustrated interaction 
respectively, while the antiferromagnetic NNN interaction always decreases NNE. Meanwhile, 
the antiferromagnetic NNN interaction can generate a larger value of NNNE in comparison to the 
case when the NNN interaction is ferromagnetic. So we can employ NNN interaction strength 
as well as dimerized interaction to realize quantum entanglement control. As for the case of 
7 7^ 1(XY model), we will present further reports in the future. 
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FIG. 1 : Diagram of the Ji — J2 — -^3 model which can be considered either as a linear chain(a) or a zigzag 
chain(b). The alternating nearest-neighbouring interaction is Ji and J2 respectively, while the next-nearest- 
neighbouring interaction is parameterized by J3 
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FIG. 2: The nearest neighbouring concurrence Ci,2 and €2,3 as a function of the reduced coupling constant 
A for various alternating values of the nearest-neighbour interaction and external magnetic field, NNN 
interaction k = 0. 
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FIG. 3: The next-nearest-neighbouring concurrence Ci^a as a function of the reduced coupUng constant A for 
various alternating values of the nearest-neighbour interaction and external magnetic field, NNN interaction 
K = 0. 
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FIG. 4: The concurrence Ci,2 and Ci^a as a function of the reduced coupling constant A for various frustrated 
values of next-nearest-neighbouring interaction, NN interaction a = (3 = 1.0. 
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